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ABSTRACT 

We  consider  in  this  paper  the  positive  0-1  polynomial  programming  (PP) 
problem  of  finding  a 0-1  n-vector  x that  maximizes  c^x  subject  to 
f(x)  < b where  c,  b ^ 0 and  f is  an  m-vector  of  polynomials  with  non- 
negative coefficients. 

Two  types  of  heuristic  methods  for  solving  PP  problems  were  developed. 

The  various  algorithms  were  tested  on  randomly  generated  problems  of  up  to 
1000  variables  and  200  constraints.  Their  performance  in  terms  of  computational 
time  and  effectiveness  was  investiaged.  The  results  were  extremely  encouraging. 
Optimal  solutions  were  consistently  obtained  by  some  of  the  heuristic  methods 
in  over  50%  of  the  problems  solved.  The  effectiveness  was  on  the  average 
better  than  99%  and  no  less  than  96.5%.  The  computational  time  using  the  heuristic 
for  PP  problems  is  on  the  average  5%  of  the  time  required  to  solve  the  problems 
to  optimality. 
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Introduction 


We  consider  In  this  paper  the  positive  0-1  polynomial  programming  (PP) 

T 

problem  of  finding  a 0-1  n-vector  x that  maximizes  c x subject  to 
f(x)  _<  b where  c,  b ^ 0 and  f Is  an  m-vector  of  polynomials  with  non- 
negative  coefficients. 

Two  main  approaches  have  been  proposed  In  the  literature  for  solving  PP 
problems.  The  first  one  Is  a linearization  method  which  converts  the  PP  prob- 
lem to  an  equivalent  linear  0-1  programming  problem  with  additional  variables 
and  constraints,  see  e.g.  [3,  4,  13].  The  second  approach  Involves  solving 
the  PP  problem  In  Its  original  form.  Methods  that  can  be  cast  Into  this  form 
are  Branch  and  Bound  [7],  Implicit  Enumeration  [10]  and  Covering  and  Generalized 
Covering  Relaxation  Algorithms  [5,  6]. 

However,  as  In  linear  Integer  programming  problems,  all  the  algorithms 
for  solving  PP  problems  suffer  from  significant  computational  limitations. 

Any  of  chose  algorithms  can  solve  only  modest  size  problems.  Motivated  by 
this  limitation  and  In  view  of  the  many  successful  heuristic  algorithms  that 
were  developed  for  linear  Integer  programs,  see  e.g.  [1,  2,  8,  9,  11,  12], 
we  will  construct  In  this  paper  some  heuristic  methods  for  PP  problems. 

The  methods  we  develop  can  be  divided  Into  two  categories.  The  first 
one  Is  a dual  approach  chat  starts  by  setting  all  variables  equal  to  one  and 
decreases  their  values,  one  at  a time,  from  one  to  zero  until  feasibility  Is 
reached  (see  also  [11,  12]  for  the  linear  case).  The  second  approach  starts 
with  a feasible  solution  to  PP  and  Improves  It  by  Increasing  Che  value  of  Che 
variables  until  no  further  Improvement  Is  possible  (see  [9]  for  the  linear 
case) . 
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The  various  algorithms  were  tested  on  hundreds  of  randomly  generated  PP 
problems  of  up  to  1000  variables  and  200  constraints.  Their  performance  In 
terms  of  computation  time  and.  In  problems  with  up  to  50  variables  and  con- 
straints, effectiveness,  l.e.,  the  percent  the  objective  function  value  of  the 
heuristic  solution  was  of  the  optimal  one,  was  Investigated.  The  computational 
results  obtained  are  extremely  encouraging.  Optimal  solutions  were  consistently 
obtained  by  some  of  the  heuristic  methods  in  over  50%  of  the  problems  solved. 
Further,  the  effectiveness  was  on  the  average  better  than  99%  and  no  less  than 
96.5%.  The  computational  time  using  the  heuristic  for  PP  problems  is  on  the 
average  5%  of  that  using  the  covering  relaxation  approach  described  In  [6]. 

2.  Preliminary  Definitions  and  Notations 

Consider  again  the  positive  0-1  polynomial  programming  PP  problem 

i 

I 

I n 

Maximize  E c x 
J-1  ^ J 


Subject  to  f^(x)  £ l~l,...,m 


Xj  ^ j-l,...,n 


where  fj(x)  are  polynomials  of  the  form 


f (x)  -la  n X 
k-l  J 
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with  N^,  any  subset  of  N ••  n}.  We  will  assume  that  the  c.  and 

Ik  j 

are  all  nonnegative  and  that  the  are  positive. 

If  we  further  assume,  without  loss  of  generality,  that  b^  > 0 (i“l,...,m) 

(since  if  not,  x = 0 is  the  only  feasible  solution),  then  we  can  normalize 
each  polynomial  function  by  dividing  all  of  its  coefficients  by  b^.  Thus 
we  can  assume  without  loss  of  generality  that  b 1.  Call  this  the  normalized 
problem. 

Now,  for  any  given  variable  , f^(x)  can  be  written  as 

f^(x)  » x^g|(x^)  + h|(x^) 

where  x'^  = (x,  ,...,x.  , ,x. x ).  We  will  refer  to  g^(x^)  as  the  deri- 

1 J-1  j+i  n 1 

vatlve  of  f^(x)  with  respect  to  the  variable  x^ . Observe  that  if  a variable 
Xj  is  Increased  from  zero  to  one  then  g|(x^)  is  the  change  in  the  left  hand 
side  of  the  i^^  constraint.  The  g|(0)  can  be  considered  as  the  consumption 
of  the  right  hand  side  incurred  by  such  an  Increase.  We  will  refer  to  g|(0) 
as  the  linearity  of  x^  in  f^(x)  and  denote  it  by  L^. 

We  will  now  briefly  motivate  the  ideas  underlying  our  heuristic  methods. 
When  solving  a PP  problem  we  would  like  to  first  set  a variable  Xj  to  one 
if  its  contribution  to  the  objective  function  is  as  large  as  possible,  while 
its  consumption  of  the  right  hand  side  is  as  small  as  possible.  Thus  we  would 


choose  to  set  x to  one  if  k*J  maximized 

J 
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Further  since  when  Is  Increased  from  zero  to  one  the  change  In  the 

left  hand  side  of  the  1^*'  constraint  Is  we  would  like  g|(x^) 

to  Involve  terms  with  small  coefficients  and  each  having  a large  number  of 
variables.  The  main  reason  for  that  Is  to  ensure  that  a further  Increase  of 
other  variables  from  zero  to  one,  in  particular  those  which  appear  In  g^(x^), 
will  consume  only  a small  amount  of  the  right  hand  side.  This  leads  to  the 
definition  of  the  weighted  linearity  of  x^  in  f^(x),  denoted  by  Wl|^,  and 
given  by 


(2) 


kel^ 


where  is  the  cardinality  of  and  g|(x^)  Is  given  by 


(3) 


g{(x3)  . E 

r<J 


kel^  '^^^Ik 


where  l|  and  are  the  appropriate  sets  defining  . 

If  on  the  other  hand  we  would  like  to  Initially  set  all  variables  equal 
to  one,  the  Infeaslblllty  of  the  1^^  constraint  will  be  given  by  (fj^(l)  -1)^  . 
If  Xj  Is  then  decreased  from  one  to  zero,  the  change  In  the  left  side  will 
be  given  by  g|(l)  and  the  Infeaslblllty  of  the  l^**  constraint  will  drop 
to  (h|(l)  - 1)'*’  where  f|(x)  Xjg|(x^)  + h|(x-^)  . Now  starting  with  x ■ 1 
and  attempting  to  solve  the  PP  problem  we  would  like  to  decrease  first  that 
variable  x^  whose  contribution  to  the  objective  function  Is  as  small  as  possible 
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while  its  consumption  of  the  right  hand  side  Is  as  large  as  possible  In  the 
violated  contralnts.  Thus  we  would  decrease  the  value  of  the  variable 
for  which  , given  by 


j 


(4) 


:i 


E (f,(i)  - i)"^  g^d) 

1=1  ^ 


Is  minimum. 

Moreover  if  g|(x^)  contains  terms  with  large  coefficients  and  which 
involve  a large  number  of  variables,  those  terms  will  be  discarded  when 
is  set  to  zero.  Thus  we  can  modify  (4)  and  choose  to  decrease  to  zero  that 
variable  for  which  WS^  , defined  by 


ws-" 


'\kl  "ikl 


1) 


kei:i 


Is  minimized. 


The  quantities  L;J  , WL"!  , and  WS'^  are  used  In  the  next  section 


j 


rJ 

d ’ 


where  the  heuristic  algorithms  for  solving  PP  are  described. 

3.  The  Heuristic  Algorithms 

We  will  present  In  this  section  two  types  of  heuristic  algorithms  for 
solving  the  PP  problem.  The  first  type  will  start  with  a feasible  solution 


5 


X 0 and  Increase  the  values  of  the  variables  one  at  a time  until  no  further 
increase  is  possible.  The  order  in  which  the  variables  are  Increased  is 
determined  using  once  their  sum  of  linearities  and  then  the  sum  of  their 
weighted  linearities.  The  second  type  of  algorithm  will  start  with  an  in- 
feasible solution  X • 1 and  use  the  Infeaslbillty  as  well  as  the  change  in 
Infeaslblllty  to  determine  variables  to  be  dropped  to  zero.  This  is  done  until 
a feasible  solution  is  reached.  Once  a feasible  solution  is  at  hand,  an  improve- 
ment procedure  will  be  applied  to  determine  whether  some  of  the  variables  dropped 
to  zero  can  be  Increased  back  to  one.  This  procedure  seems  to  play  a crucial 
role  in  the  excellent  performance  of  the  second  type  of  algorithms  on  the  PP 
problem. 

The  six  different  heuristic  algorithms  can  now  be  formally  stated  as  follows 
Algorithm  1; 

Step  0;  Start  by  setting  x » 0 , Iq  “ {l,...,n}  , 

I “ 0 and  M » {l,...,m}  . 

r 


Step  1; 


Let 


j be  an  index  k e maximizing 


Step  2;  Check  whether  by  increasing  from  zero  to  one  any  of  the 

constraints  will  be  violated.  If  yes  go  to  step  3;  otherwise 
go  to  step  4. 

Step  3;  Set  1^  = 1^  U { j } and  1^  » ' 


If  Iq  = (9  go  to  step  5;  otherwise  go  to  step  1. 


Step  4: 


Set  X 


■j 


1 , subtract  from  both  sides  of  the  i^^  con- 


straint, and  normalize  the  problem,  thereby  determining  modified 
J c-t-  T = T VI  fj  T « 0 . If  I = 0 , 


a^j  and 


Set  Iq  = Iq  U (Ip\{j})  , Ip 


go  to  step  5;  otherwise  if  the  k 
set  M “ M\{k}  and  go  to  step  1 


th 


constraint  is  redundant 


Step  5; 


Terminate  with  x. 


Algorithm  I can  be  modified,  using  the  weighted  linearities,  to  produce 
Algorithm  II; 

Same  as  algorithm  I except  replace  step  1 by 


Step  1;  Let  j be  an  index  k e Iq\Ip  maximizing 


Z wl' 

ieM 


k • 
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The  next  four  algorithms  will  start  wx.h  the  optimal  solution  x » 1 to 
the  relaxation  of  the  PP  problem  In  which  the  polynomial  constraints  are  dropped 
and  proceed  towards  feasibility. 

Algorithm  III 

Step  0;  Start  with  x = 1 . If  x Is  feasible  for  PP,  terminate.  Otherwise 
set  J » {l,...,n}  . 


Step  1:  Determine  the  set  I of  indices  1 for  which  f^(x)  > 1 


Go  to  step  2. 


Step  2;  Let  k be  an  Index  j in  J minimizing 


'± 


E (f/x)  - l)"^  gJCx^) 
lei 


Set  x^  = 0 . If  X is  feasible,  terminate  with  x;  otherwise 
replace  J by  J\{k}  and  go  to  step  1. 

Step  2 of  algorithm  III  can  be  further  modified  resulting  with  algorithm 
IV  as  follows: 

Algorithm  IV 

Same  as  algorithm  III  except  that  in  step  2 we  let  k be  an  Index  J 
that  minimizes 


WS-" 


c . 


I (f,(x) 
iei 


- 1) 


ir’  ir 


rei-^ 


where  N,  and  are  defined  in  (3). 

ir  i 

Observe  that  in  algorithms  III  and  IV  whenever  a variable  is  set 

to  zero  the  algorithms  return  to  step  1 and  calculate  the  next  best  candidate 
to  be  set  to  zero.  The  following  two  algorithms  are  similar  to  algorithms 
III  and  IV,  however  the  variables  are  ordered  once  and  then  set  to  zero  one 
at  a time  until  feasibility  is  reached. 

Algorithm  V: 

Step  0;  Start  with  x = 1 . If  x is  feasible  for  PP,  terminate.  Otherwise 
go  to  step  1. 


Step  1: 


Let  j 

X n 

Ji  < if  s" 
equal  to  zero. 
Terminate  with 


be  a permutation  of  the  first  n integers  such  that 

■I  jif  j-{  ji, 

^ < S or  if  i < k and  S ^ - S . Set  x^ 
in  order,  one  at  a time,  until  feasibility  is  reached. 

X . 


Algorithm  VI; 

Same  as  algorithm  V except  the  variables  are  ordered  in  Increasing  order 
of  WS^  . 

The  solution  x determined  by  any  of  the  algorithms  III  to  VI  can  be 
further  Improved.  This  Improvement  was  found  to  play  a significant  role  in 
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the  effectiveness  of  the  methods  of  the  second  type.  To  motivate  this  pro- 
cedure observe  that  a decision  to  drop  the  value  of  a given  variable  from  one 


to  zero  may  very  well  be  reversed  after  other  variables  having  the  value  one 
are  forced  to  zero.  Thus  after  determining  a solution  x , we  would  like 
to  check  whether  any  of  the  variables  dropped  to  zero  can  be  increased  back 
to  one  without  violating  feasiblity.  This  can  be  easily  done  using  the  follow- 
ing procedure: 

Modification  Procedure: 


Step  0: 


Start  with  x — the  solution  obtained  by  any  of  the  heuristic 
algorithms  III  to  VI.  Let  PP  be  the  problem  obtained  from  PP  after 
substituting  the  values  of  all  variables  x^  that  are  equal  to  one 
in  X . (Thus  PP  involves  only  those  variables  whose  value  is  zero 
in  X. ) 


Step  1:  Use  heuristic  method  II  for  W,  terminating  with  x.  Decrease  to  one 

the  values  of  those  variables  in  the  x from  step  C that  are  equal 

A* 

to  one  in  x , and  terminate  with  this  new  x . 


I 

i 


The  above  improvement  process  can  be  easily  modified  to  incorporate  cases 
where  we  would  like  to  test  whether  a decrease  of  another  variable  from  one 
to  zero  can  Improve  the  heuristic  solution.  This  is  done  by  eliminating  the 
substitution  of  that  variable  in  PP  in  step  0 of  the  modification  process. 


r 

I 


> 4 . Computational  Results 

The  six  algorithms  were  coded  in  Fortran  IV  and  Implemented  on  an  IBM 
370/168.  The  performance  of  the  algorithms  was  tested  on  a large  number  of 
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randomly  generated  problems  and  was  measured  in  terms  of  computational  time 
and  effectiveness,  l.e.,  the  percent  the  objective  function  value  of  the  heuris- 
tic solution  was  of  the  optimal  one.  For  the  larger  problems  with  up  to  1000 


variables  and  200  constraints,  for  which  an  optimal  solution  was  not  sought, 

the  algorithms  were  compared  according  to  their  relative  effectiveness. 

We  have  tested  the  performance  of  the  heuristic  algorithms  in  terms  of 

five  parameters,  viz.,  the  number  n of  variables,  the  number  m of  constraints, 

the  maximum  number  k of  terms  per  constraint,  the  maximum  number  v of 

variables  in  each  term,  and  the  degree  of  tightness  a of  the  constraints. 

The  data  for  each  PP  problem  was  randomly  generated  as  follows.  The 

cost  vector  c was  determined  by  setting  c^  “ 0 and  ^j+2.  * ^ ^ where 

p is  randomly  chosen  from  [0,10].  The  coefficients  a are  randomly  chosen 

^ m 

between  0 and  10,  while  the  right  hand  side  b was  set  equal  to  a • E a , 

j=l 

with  a equal  to  0.3,  0.5  or  0.9.  The  number  of  terms  in  each  constraint 
was  randomly  chosen  between  1 and  k and  the  number  of  variables  in  each 
term  was  randomly  chosen  between  1 and  v. 

All  possible  combinations  obtained  by  varying  n and  m among  30, 

40,  50  and  a among  0.3,  0.5  and  0.9  were  tested.  From  the  results  obtained 
for  the  270  problems  that  were  run,  ten  for  each  combination  of  n,  m and 
a,  we  can  conclude  that  method  II  dominates  method  I,  and  methods  III  and  IV 
are  uniformly  superior  to  methods  V and  VI  both  in  terms  of  efficiency  and 
in  computational  time.  There  was  no  clear-cut  choice  between  methods  III  and 
IV.  All  problems  were  initially  solved  with  the  original  algorithms  without 
adding  the  modification  procedure.  The  results  obtained  by  method  II  were 
extremely  good  and  the  effectiveness  was  never  below  96. 5X.  The  performance 


of  methods  111  and  IV  was  not  as  good  and  In  some  cases,  especially  when  Che 
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constraints  were  tight,  the  effectiveness  sometimes  fell  to  94-95%.  However, 
both  methods  III  and  IV  required  less  computation  time  to  produce  their  heuristic 
solutions.  After  employing  the  modification  to  methods  III  and  IV  the  solutions 
were  uniformly  improved  and  many  of  them  reached  the  optimal  value.  In  most 
cases,  except  when  the  constraints  were  very  loose,  the  total  computational 
time  required  by  the  modified  methods  III  and  IV  exceeded  that  of  method  II. 

Of  the  many  randomly-generated  problems  solved,  method  II  reached  optimality 
in  about  50%  of  the  problems,  while  methods  III  and  IV  did  so  in  about  30% 
of  the  problems.  After  employing  the  modification  process  on  methods  III  and 
IV,  they  reached  optimality  in  about  60%  of  the  problems  solved. 

The  influence  of  the  number  of  variables,  number  of  constraints  and  tight- 
ness of  the  constraints  on  computational  time  and  effectiveness  in  algorithm  II 
and  the  modified  algorithms  III  and  IV  is  summarized  in  Table  1 and  Figure  1. 

The  slopes  a and  the  intercepts  b for  each  line  in  Figure  1 were  obtained 
using  averages  from  all  runs  in  which  one  of  the  factors  was  kept  constant  and 
all  others  were  changed  in  their  given  ranges. 

The  effectiveness  of  all  three  methods  was  not  affected  by  either  the 
number  of  constraints  or  the  number  of  variables.  However,  as  shown  in  Figure  1 
the  effectiveness  decreases  for  method  II  and  increases  for  methods  III  and  IV 
with  a decrease  in  constraint  tightness.  Computation  time,  however.  Increases 
for  all  three  methods  with  an  Increase  in  the  number  of  variablea  and/or  con- 
straints. This  increase  in  computation  time  is  strongly  Influenced  by  the 
tightness  of  the  constraints.  Indeed  method  II  is  faster  for  tight  constraints, 
while  methods  III  and  IV  are  faster  for  loose  constraints.  By  comparing 

j 

the  performance  of  the  three  methods  we  observe  that  for  a » 0.3  and  0.5,  I 
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The  computational  time  reported  is  the  true  execution  time  (total  CPU  time  minus  loading  input/output  and 
overhead  time)  given  in  milliseconds  i.e.  l/lOOO  of  a second. 


method  II  Is  almost  uniformly  faster  than  methods  III  and  IV.  However,  for 
a “ 0.9  methods  III  and  IV  become  faster  than  method  II,  especially  for  the 
larger  problems.  In  fact,  for  m » 0.9  methods  111  and  IV  without  the  modi- 
fication process  are  as  effective  as  method  II.  Moreover,  they  require  sig- 
nificantly less  computation  time  to  produce  the  heuristic  solutions  than  Is 
required  by  method  11.  These  results  are  summarized  in  Table  2.  Indeed  in 
Table  2 we  observe  that  the  effectiveness  of  the  three  methods  is  about  the 
same  and  that  methods  III  and  IV  are  almost  three  times  faster  than  method  II. 

In  Table  3 and  Figure  3 we  exhibit  the  Influence  of  a change  in  the  maxi- 
mum number  of  terms  k In  each  constraint  and  the  maximum  number  of  variables 

V in  each  term  on  both  effectiveness  and  computation  time. 

The  computation  time  Is  Increased  for  all  three  algorithms  with  an  Increase 
in  either  k or  v , however  with  a much  larger  slope  for  method  II  than  for 
methods  III  and  IV.  Again  the  effectiveness  of  neither  of  these  methods  is 
superior  to  the  other  two.  It  is  interesting  to  note  that  when  either  k or 

V is  very  small  (average  of  2)  the  effectiveness  of  all  three  methods  is 
worst. 

The  performance  of  methods  II  and  modified  methods  III  and  IV  was  then 
tested  on  some  large  problems  of  up  to  1000  variables  and  200  constraints  that 
were  generated  in  a similar  way  to  those  generated  before.  No  attempt  has 
been  made  to  find  optimal  solutions  for  those  problems  because  of  the  excessive 
computation  time  It  would  require.  The  performance  of  the  three  methods  was 
compared  by  means  of  their  effectiveness  relative  to  the  best  heuristic  obtained 
and  by  time  relative  to  the  worst  time  obtained.  The  results  were  averaged 
from  ten  runs  and  are  summarized  in  Table  4.  i 
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ZAT  equals  the  percentage  the  average  time  required  is  of  the  worst  time. 


From  Table  4 we  observe  that  the  relative  effectiveness  of  modified 
methods  HI  and  IV  seem  to  improve  slightly  with  an  increase  in  the  number 
of  variables.  It  is  Interesting  to  note  that  the  relative  effectiveness  of  the 
three  methods  remains  nearly  constant,  ranging  between  98.39  to  99.99%.  The 
relative  time  on  the  other  hand  decreases  almost  linearily  for  modified  methods 
III  and  IV.  This  makes  them  more  attractive  for  large  size  problems. 


I 

i 

1 


21 


ACKNOWLEDGEMENT 


The  author  wishes  to  acknowledge  the  programming  assistance  of 
Mr.  W.  Vaessen.  The  author  is  also  greatful  to  Professor  Daniel  Granot 
for  valuable  discussions  on  the  subject,  and  to  Professor  A.  Veinott  for 
his  many  helpful  remarks  on  an  earlier  draft  of  the  paper. 


References 


[1]  BAl.AS,  E.  and  C.  H.  MARTIN,  "Pivot  and  Cornplcment  - A Heuristic  for 

0-1  Programming",  in  MSR  Report  No.  ^14,  Carnegie  Mellon  University, 
February  1978. 

[2]  FAALAND,  U.  H.  and  F.  S.  HIELiER,  "Interior  Path  Methods  for  Heuristic 

Integer  Programming  Procedures",  Technical  Report  1173,  Department  of 
Operations  Research,  Stanford  University,  February  1977. 

[3]  GLOVER,  F.  and  E.  WOOLSEY,  "Further  Reduction  of  Zero-One  Polynomial 

Programming  Problems  to  Zero-One  Linear  Programming  Problems", 
Operations  Research,  Vol.  21,  No.  1 (1973),  pp.  141-161. 

[4]  GLOVER,  F.  and  E.  WOOLSEY,  "Converting  the  0-1  Polynomial  Programming 

Problem  to  a 0-1  Linear  Program",  Operations  Research,  Vol.  22,  No.  1 
(1974),  pp.  180-182, 


[5]  GRANOT,  D.  and  F.  GRANOT,  "Generalized  Covering  Relaxation  for  0-1 

Programs",  Technical  Report  SOL  78-16,  Stanford  University,  June  1978. 

[6]  GRANOT,  D.,  F.  GRANOT  and  J.  KOLT.BERG,  "Covering  Relaxation  in  Monotone 

0-1  Programming",  Submitted  to  Management  Science. 

[7]  HAMMER,  P.  L. , "A  B-B-B  Method  for  Linear  and  Nonlinear  Bivalent 

Programming",  Operations  Research,  Statistics  and  Economics  Mimeo- 
graph Series  No,  48,  Technlon,  May  1969. 

[8]  HILLIER,  F,  S.,  "Efficient  Heuristic  Procedures  for  Integer  Linear 

Programming  with  an  Interior",  Operations  Research,  17  (1969), 
pp.  600-637. 

[9]  KOCKENBERG,  G.  A.,  B.  A.  MCCOIL  and  F.  P.  WYMAN,  "A  Heuristic  for 

General  Integer  Programming",  Decision  Sciences,  Vol.  5 (1974), 
pp.  36-44, 

110]  LAUGHHUM,  D.  L. , "Quadratic  Binary  Programming  with  Application  to 
Capital  Budgeting  Problems",  Operations  Research,  Vol.  18,  No.  3 
(1970),  pp.  454-461. 


23 


[11]  SENJrU,  S.  and  Y.  TOYODA,  "An  Approach  to  Linear  Progrannnlng  with  0-1 
Variables",  Management  Science,  15  (1968),  pp.  B196-207 . 


[12]  TOYODA,  Y.,  "A  Simplified  Algorithm  for  Obtaining  Approximate  Solutions 

to  0-1  Programming  Problems",  Management  Science.  21  (1975),  pp.  1417-1427. 


[13]  WATTERS,  L.  J.,  "Reduction  of  Integer  Polynomial  Programming  Problems 

to  Zero-One  Linear  Programming  Problems",  Operations  Research,  Vol.  15 
(1967),  pp.  1171-1174. 


[14]  ZANAKIS,  S.  M.,  "Heuristic  0-1  Linear  Programming:  An  Experimental 

Comparison  of  Three  Methods",  Management  Science,  24  (1977),  pp.  91-104. 


24 


UNCLASSIFItD 

SCCUHITV  CLAtSiriCATION  OP  THIS  PAOC  fWAan  Oats  Cntaf^ 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLBTINa  FORM 

I.  RCClRlENT‘1  catalog  NUMBER 

4-  TiTue  (md  Subf/rl«> 

EFFICIENT  HEURISTIC  ALGORITHMS  FOR  POSITIVE 

0-1  POLYNOMIAL  PROGRAMMING  PROBLEMS 

S.  TYRE  or  REPORT  A PERIOO  COVERED 

TECHNICAL  REPORT 

t.  PERFORMING  ORO.  REPORT  NUMBER 

7.  AUTHONfal 

Frieda  Granot 

t.  CONTRACT  OR  GRANT  NUMBERfaf 

N00014-76-C-0418‘^ 

• . PERPONMIMG  ONOANIZATION  NAME  AND  ADDRCft  / 

DEPARTMENT  OF  OPERATIONS  RESEARCH  > 

STANFORD  UNIVERSITY,  STANFORD,  CALIF. 

H.  CONTROLLING  OPPICE  NAME  ANO  AOORCU 

OPERATIONS  RESEARCH  PROGRAM  CODE  434 

OFFICE  OF  NAVAL  RESEARCH 

ARLINGTON,  VIRGINIA  22217 

U.  REPORT  DATE 

Auaust  1978 

IS.  number  of  pages 

24 

l«.  MONITORING  AGENCY  NAME  A AOOREIVIf  MIHtmt  hum  CMiCralfRiE  OHItu) 

IS.  security  CLASS,  faf  Mia  raparff 

UNCLASSIFIED 

If.  OISTRItUTfOM  STATCMCNT  (dt  thU  Kdpdft) 

APPROVED  FOR  PUBLIC  RELEASE  AND  SALE;  DISTRIBUTION  IS  UNLIMITED. 

17.  DISTRIBUTION  STATEMENT  (ut  Utu  abairaci  anlaraV  hi  Bluuk  M,  II  EllfaaanI  hum  Rapafi; 

20  on  Contract 

This  report  also  issued  as  Technical  Report  #78- 
Contract  F 44620-74-C-0079,  Air  Force. 

IS.  KEY  WORDS  fCatillnua  an  rararaa  alEa  If  naaaaaaqr  ant  lEanHIp  Sr  SlaaS  tumhui) 

Integer  Polynomial  Programming  Algorithms 

Integer  Programming  Polynomial  Prograinning 

Heuristic  Algorithms 

to.  abstract  (CuMImtu  an  rararaa  alOa  If  naanaamr  anO  UrnuHB  kf  kluuk  uumktk) 

SEE  REVERSE  SIDE 

DO  I jAiTn  1473  totrioN  OP  I NOV  m it  omoliti 

t/H  titt-ei«>t«ot  I 


MCWMTV  ^AM 


vsa 


UNCLASSIFIED 

iKCUWiTv  ccawficatiom  of  Thu  nw>—  Om* 

EFFICIENT  HEURISTIC  ALGORITHMS  FOR  POSITIVE  0-1  POLYNOMIAL  PROGRAMMING 
PROBLEMS 
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We  consider  in  this  paper  the  positive  0-1  polynomial  programming  (PP) 
problem  of  finding  a 0-1  n-vector  x that  maximizes  c^x  subject  to 
f(x)  £ b where  c,  b £ 0 and  f is  an  m- vector  of  polynomials  with 
non-negative  coefficients. 

Two  types  of  heuristic  methods  for  solving  ^(PP)  problems  were  developed. 

The  various  algorithms  were  tested  on  randomly  generated  problems  of  up  to 
1000  variables  and  200  constraints.  Their  performance  in  terms  of  compu- 
tational time  and  effectiveness  was  investiaged.  The  results  were  extremely 
encouraging.  Optimal  solutions  were  consistently  obtained  by  some  of  the 
heuristic  methods  in  over  50%  of  the  problems  solved.  The  effectiveness 
was  on  the  average  better  than  99%  and  no  less  than  96.5%.  The  computa- 
tional time  using  the  heuristic  for  PP  problems  is  on  the  average  5% 
of  the  time  required  to  solve  the  problems  to  optimality. 
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